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Addition  of  Vertical  Velocity  to  a  One-Dimensional 
Aerosol  and  Trace  Gas  Model 


I.  INTRODUCTION 

The  NRL  aerosol  model  includes  all  aerosol  processes  believed  to  be  important  in  the 
Marine  Boundary  Layer  (MBL).  These  processes  can  be  divided  into  three  classes  with  respect 
to  their  dependence  on  meteorology:  (1)  those  processes  that  are  only  incidentally  dependent  on 
meteorology,  (2)  those  that  are  driven  almost  entirely  by  the  meteorology  and  (3)  those  for  which 
both  the  meteorology  and  the  aerosol  properties  are  crucial.  In  the  first  class  are  processes  such 
as  coagulation,  dififusional  growth  of  core  material  by  condensation  of  gas-phase  reaction 
products  (such  as  sulfate),  gravitation  settling,  and  nucleation  of  new  particles.  These  processes 
are  determined  at  the  molecular  mean-free-path  and  aerosol  size  scales.  The  ambient  relative 
humidity  affects  these  processes  indirectly  through  the  change  in  aerosol  size  by  uptake  of  water 
(for  hygroscopic  aerosol  components).  The  strength  of  atmospheric  turbulence  can  have  a  small 
inertial  effect  on  coagulation.  In  the  second  class  of  processes,  such  as,  horizontal  advection, 
vertical  mixing  of  aerosols,  and  exchange  with  the  free  troposphere,  meteorological  conditions 
play  the  dominant  role.  Other  processes  are  crucially  dependent  on  both  the  ciurent  meteorology 
and  microphysical  properties  of  the  aerosol.  Examples  of  these  are  the  various  precipitation 
scavenging  mechanisms,  addition  of  aerosol  mass  during  cloud  processing  (liquid-phase 
chemistry),  surface  deposition  and  generation  of  sea-salt  aerosol  at  the  sea  surface. 

It  is  clear  from  the  above  that  it  is  impossible  to  model  the  spatial  and  size  distribution  of 
aerosols  in  the  MBL  without  having  an  accurate  concurrent  description  of  the  meteorological 
fields  (wind,  turbulent  mixing,  RH,  temperature,  clouds  and  precipitation  fields).  The  most 
robust  and  accurate  aerosol  models  are  sectional  models  where  the  aerosol  size  distribution  is 
divided  into  a  number  of  size  sections  with  each  section  containing  a  number  of  (chemical) 
components.  The  NRL  sectional  model  eliminates  stiffness  and  numerical  diffusion  caused  by 
calculating  condensation  and  evaporation  of  water  vapor  during  movement  of  the  aerosol 
through  relative  humidity  gradients  and  greatly  reduces  numerical  diffusion  for  growth  of  core 
material  by  various  gas-to-partiele  conversion  mechanisms  (Fitzgerald  et  al.[1998a,b]  and 
Gelbard  et  al.  [1998]). 

The  NRL  (1-D  aerosol  model  is  typically  run  with  37  size  sections  (between  0.005  and  15 
pm  dry  radius)  and  four  components,  sulfate,  sea  salt,  an  insoluble  component  (dust),  and  water 
in  equilibrium  with  the  soluble  components  at  ambient  RH.  To  fully  integrate  an  aerosol  model 
into  a  3-D  meteorological  model  would  require  carrying  and  calculating  over  100  aerosol 
variables  at  each  grid  point  of  the  model  in  addition  to  the  meteorological  variables.  While  the 
rapid  increase  in  computing  speed  may  make  this  possible  in  the  future,  or  even  the  present, 
provided  enough  computing  time  is  available,  we  have  settled  on  a  less  ambitious  model 
integration  scheme  that  uses  the  output  of  a  high  resolution  mesoscale  model  to  drive  our  1-D 
aerosol  model. 

The  model  integration  being  pursued  is  the  following:  A  high-resolution  mesoscale 
model,  in  our  case  the  Navy’s  COAMPS  model,  is  run  over  the  geographical  region  and  time  of 
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interest  and  the  meteorological  variables  at  each  grid  point  are  stored  at  specified  time  intervals, 
typically  every  hour.  In  addition  to  the  normal  meteorological  fields,  some  internally  generated 
parameters,  such  as  the  turbulent  mixing  coefficients  generated  by  the  boundary  layer 
parameterization  used  in  the  model,  are  also  saved.  From  the  horizontal  wind  field,  a  mean  MBL 
back  trajectory  from  any  point  of  interest  is  generated.  The  NRL  1-D  aerosol  model,  here 
considered  as  a  vertical  column,  is  advected  along  this  trajectory.  The  metrologi'cabdata  required 
by  the  aerosol  model  is  extrapolated  using  the  grid  points  that  lie  closest  to  the  points  on  the 
back-trajectory.  Because  of  the  well-mixed  character  of  the  MBL,  the  concept  of  a  1-D  column 
moving  with  the  mean  wind  speed  in  the  MBL  is  a  reasonable  approximation.  The  extension  of 
the  column  into  the  lower  free  troposphere  (FT)  to  simulate  MBL-FT  exchange  is  more 
questionable  (because  of  possible  velocity  gradient  across  the  boundary  and  lack  of  mixing).  In 
some  large-scale  experiments,  an  exchange  velocity  is  derived  from  the  measured  exchange  of 
ozone  or  water  vapor;  in  which  case  we  terminate  the  column  at  the  top  of  the  MBL  and  use  the 
measured  exchange  velocity.  For  the  case  where  we  have  only  the  numerically  generated  data, 
the  exchange  coefficients  themselves  define  the  extent  of  the  MBL  and  the  exchange  with  the  FT 
is  determined  by  the  magnitude  of  these  coefficients,  plus  any  exchange  that  results  from 
synoptic-scale  vertical  velocities  generated  by  the  mesoscale  model. 

The  mesoscale  model  calculates  (sub-grid)  turbulent  mixing  by  the  boundary  layer 
parameterization  used  in  the  model.  Large-scale  subsidence/lifting  (large  compared  to  the  grid 
scale)  are  resolved  by  the  numerical  solutions  to  the  fluid  mechanical  equation.  The  turbulent 
mixing  coefficients  describe  the  statistical  effect  of  mixing  due  to  eddies  of  the  smallest  sizes  up 
to  and  including  the  so-called  “large  scale  eddies”  which  have  dimensions  comparable  to  the 
height  of  the  MBL.  The  characteristic  length  which  separates  the  largest  eddies  from  the  large- 
scale  vertical  velocities  generated  by  the  model  is,  of  course,  not  well  defined.  But  the  same 
division  of  the  two  mechanisms,  which  results  from  the  mesoscale  model  treatment,  is  used  in 
the  aerosol  calculations. 

In  the  original  aerosol  model  (Fitzgerald  et  al.,  1998a,  1998b;  Gelbard  et  al.,  1998)  large- 
scale  vertical  motions  were  not  included.  Locally,  the  effects  of  turbulent  mixing  in  the  MBL  are 
much  stronger  than  large-scale  divergence  which  can  often  be  neglected.  Over  longer  periods  of 
time  and  for  exchange  with  the  free  troposphere,  the  large-scale  vertical  motions  would  appear  to 
be  important.  If,  for  no  other  reason,  inclusion  of  large-scale  vertical  velocities  is  necessary  to 
evaluate  just  how  important  vertical  velocities  are  relative  to  turbulent  mixing  and  when,  if  ever, 
the  large  scale  motions  can  be  ignored.  Availability  of  the  output  of  the  mesoscale  model  gives  a 
set  of  vertical  velocities  and  mixing  coefficients  that  should  be  consistent  and  realistic.  The 
possible  importance  of  vertical  velocity  is  suggested  by  the  observations  (J.S.  Reid,  private 
communication),  which  found  sea-salt  aerosol  in  the  outflow  regions  from  tropical  disturbances. 
Whether  or  not  the  updrafts  that  are  responsible  for  these  events  are  on  a  large  enough  horizontal 
scale  to  be  captured  by  the  model  has  yet  to  be  evaluated. 

This  report  describes  how  vertical  motions  are  included  in  the  current  1-D  aerosol  model. 
The  task  of  including  vertical  velocities  was  more  difficult  than  originally  anticipated,  largely 
because  the  amount  of  mass  added  or  subtracted  to  a  vertical  cell  (due  to  horizontal 
convergence/divergence)  must  exactly  balance  the  difference  in  that  which  flows  through  the 
upper  and  lower  cell  boundaries  for  each  size  section  in  the  case  when  there  is  no  vertical 
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gradient  in  the  mass  densities  of  the  sections.  The  condition  that  there  he  no  change  in  the  mass 
in  a  vertical  cell  when  there  was  no  vertical  mass  gradient  for  all  convergent/divergent  cases  of 
air  flow  is  the  primary  test  for  judging  whether  or  not  the  algorithm  is  correct.  Any  imbalance 
gives  a  spurious  source/sink  of  particles  that  can  causes  a  large  excess/deficit  of  aerosol  mass 
with  time.  Of  course,  if  there  is  any  gradient  in  the  mass  density,  that  gradient  must  be  moved 
up/down  depending  on  the  vertical  velocity. 

II.  ADDITION  OF  VERTICAL  VELOCITY  TO  A  ONE-DIMENSIONAL  MODEL 


A.  Basic  Scheme 

Aerosol  mass  in  each  size  section  is  conserved  during  displacements  resulting  from  air 
motion  and  gravitational  settling.  If  the  aerosol  mass  density  of  particles  in  size  section,  i,  is  mj , 
then  conservation  of  mass  requires 

^  +  V.(m,v,)=0  (1) 

where  is  the  sum  of  the  fluid  and  settling  velocity  wf 


V,  =v^+wf  (2) 

We  will  now  drop  the  size  index,  i,  with  the  understanding  that  the  mass  in  each  size  section  is 
conserved  individually  during  displacements.  Integration  of  Eq.  (1)  over  a  vertical  column  of 
volume  V  with  base  of  area  A  gives  the  change  of  the  total  mass  M  in  the  column 

•  {fniv^V  =  -  •  dS  (3) 

dt  y 


If  the  column  spans  one  vertical  cell  then  M  becomes  the  mass  of  one  size  section  in  that 
cell  and  Eq.  (3)  becomes 


dM 


dt 


—  =  -w,m,A  +  w^m^A  -  ^ 

hor 


mv^  •  dS 


(4) 


where  w,  and  are  the  vertical  velocity  components  at  the  top  and  bottom  of  the  cell  and 

consist  of  both  a  fluid  and  gravitational  settling  component,  the  latter  is  different  for  each  size 
section;  i.e., 

w,  =  +  wf  (where  i  here  refers  to  size  section  i). 
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The  horizontal  surface  integral  (Eq.  (4))  represents  the  amount  of  mass  that  enters  the  cell 
through  the  sides  of  the  column.  For  a  1-D  model,  m  is  the  same  at  all  points  around  the 
periphery  of  the  column  so  that  Eq.  (4)  can  be  written  as 

=  +  •dS  (5) 

her 

Where  (m)  is  the  vertically  averaged  mass  density  for  the  section. 

The  above  conceptual  picture  cannot  be  precisely  correct  for  a  truly  1-D  model.  In  a  1-D 
model  there  can  be  no  flow  into  the  sides  of  the  column,  and  neither  can  there  be  a  divergence  in 
the  horizontal  velocity;  i.e.,  Wt  cannot  be  different  than  Wb.  To  those  who  object  to  the  above 
conceptual  picture  they  should  view,  what  is  called,  the  horizontal  surface  integral  in  Eq.  (5)  as  a 
source  term  that  must  be  added  to  the  1-D  layer  to  balance  any  divergence  in  the  vertical 
velocity,  imposed  externally  on  the  layer.  In  practice  we  will  impose  a  vertical  velocity  gradient 
calculated  by  the  3-D  meteorological  model  at  the  location  of  the  column.  This  source  term 
compensates  for  the  convergence/divergence  in  the  fluid  flow  by  bringing  in  particles  from  the 
same  layer  (in  a  manner  which  consistent  with  the  1-D  hypothesis). 


Mass  conservation  for  an  incompressible  fluid  is  given  by 

V*v'^=0  (6) 

Integrating  Eq.  (6)  over  the  volume  of  the  column  (within  a  cell)  as  done  above  for  aerosol  mass 
yields 


-w^A  +  w^A-  •  dS  =  0  (7) 

hor 

Equations  (7)  and  (5)  can  be  combine  to  eliminate  all  references  to  the  horizontal  flow  and  area 
of  the  side  of  the  column 

J_  ^^ceii  _  +  wf  )/n,  -  (w/  +  wf  )m4  ]  +  (m)(w/  -  w/  )  (8) 

A  dt 

This  equation  satisfies  the  necessary  condition  that  if  there  is  no  gradient  in  the  mass,  m^=m^  = 
(m) ,  then  there  is  no  change  in  the  mass  in  the  cell,  even  if  there  is  a  divergence  in  the  vertical 

velocity;  provided  of  course,  that  the  settling  velocity  is  the  same  across  the  cell  wf  =  .  This 

latter  condition  will  hold  only  in  the  case  that  the  RH  is  constant  with  height.  If  one  objects  to 
the  mixture  of  1-D  and  3-D  models  used  in  the  above  derivation,  Eq.  (8)  can  be  justified  alone  on 
the  bases  of  the  following  necessary  condition:  in  the  1-D  model  there  can  be  no  change  in  the 
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mass  in  the  cell  if  there  is  no  gradient  in  the  mass  density  even  if  there  is  a  gradient  in  the 
vertical  velocity  (in  the  case  of  no  RH  gradients  and  incompressible  flow). 

The  above  scheme  is  implemented  numerically  by  assigning  a  fluid  velocity  to  cell,  n, 
and  testing  the  net  velocity  (wf  +  wf  *) ,  for  a  given  section  denoted  by  i.  Since  the  velocity 

changes  across  the  vertical  cell,  and  can  even  change  signs  within  the  cell,  assigning  a  single 
velocity  to  a  cell  introduces  certain  (finite  difference)  errors  which  can  be  severe  if  they  are  not 
handled  correctly  in  the  logical  structure  of  the  program.  Some  of  these  complications  will  be 
alluded  to  later.  For  net  upward  flow  in  the  n*'’  cell  and  both  adjacent  cells,  the  change  in  mass 
of  the  n*''  cell  is 


motion  in  adjacent  cells 


+  <  ]+  (W/  -  )  (11) 

The  logic  is  actually  more  complicated  because  adjacent  cells  can  have  velocities  in  different 
direction  so  that  there  are  four  cases  when  the  net  velocity  is  up  in  cell  n  and  four  cases  when  the 
net  velocity  is  down  in  cell  n,  giving  a  total  of  eight  cases  (See  Figure  1).  For  the  first  terms  of 
the  above  two  equations,  this  logic  is  obvious  and  straight  forward,  since  it  is  a  matter  of  looking 
at  the  net  flow  cell  by  cell  and  removing  the  correct  amount  of  mass  and  adding  it  to  the  correct 
adjacent  cell.  It  is  less  obvious  but  crucial  that  the  second  term  is  calculated  such  that  the  airflow 
into  or  out  of  the  given  cell  is  consistent  with  the  boundaries  across  which  we  move  the  particles 
using  the  net  particle  velocity.  If  this  is  not  done  properly,  a  small  imbalance  can  act  as  an 
erroneous,  and  potentially,  large  source  (or  sink)  of  mass  that  accumulates  with  time.  Figure  1 
illustrates  how  the  net  fluid  flow  (expression  in  parenthesis  in  the  last  term)  is  calculated  for  the 
eight  cases.  Since  the  settling  velocity  is  a  function  of  size,  this  procedure  is  done  separately  for 
each  section.  For  gases  only  the  air  velocities  are  required,  i.e.,  w®  =  0 .  For  particles  the  net 
velocity,  air  velocity  plus  gravitational  settling,  determines  the  direction  and  boundary  which 
must  be  considered,  e.g.,  one  section  may  be  moving  up  and  another  downward. 

We  have  assumed  here  that  the  atmosphere  is  an  incompressible  fluid;  i.e.,  no  density 
variations.  Appendix  A  shows  that  the  vertical  gradient  of  density  in  the  atmosphere  can  be 
taken  into  account.  The  interpretation  of  the  results,  and  hence  testing  of  the  algorithms  to  verify 
the  accuracy  of  the  algorithms,  is  more  complicated  for  the  compressible  case.  Therefore  we 
have  here  assumed  that  the  atmosphere  is  an  incompressible  fluid. 


1  dM 


A  dt 


Likewise  for  downward 


1  dM. 


A  dt 


- = -[(«'/ 
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w(rH-1)<0 


i 


w(nf1)>0 


w(n)>0 


w(n-1)>0 


w(rH-1)<0 


w(n)>0 


w(r>1)>0 


4 

f 


w(r><-1)>0 


1  w(n)>0 
w(i>1)<0 


w(n)>0 


w(i>1)<0 


A(Wa)=Wa(n+l)-Wa{n)  A(Wa)=-Wa(n)  A(Wa)=-Wa(n)  A(Wa)=-Wa(n)-a(n-l) 

+Wa(n+1)-Wa(n-1) 


Figure  1.  If  w  is  the  net  vertical  velocity  including  both  gravitational  lettling  and 
vertical  air  velocity  and  Wa  is  the  air  velocity,  then  the  added  air  (velocity)  consistent 
with  the  vertical  motion  of  the  particle  in  the  given  cell  and  section  is  given  by  A(Wa). 
For  downward  motion  of  the  particle  in  cell  n,  A(Wa)  is  the  appropriate  airflow  into  or 
out  of  the  cell  n  for  a  given  section,  where  positive  flow  is  downward.  There  is  a 
similar  diagram  when  the  flow  in  cell  n  is  upward. 


B.  Aerosol flux  in  terms  of  density  functions. 

In  the  analysis  above  leading  to  Eq.  (8),  the  particles  in  a  section  are  assumed  to  be 
monodisperse  whereas  in  the  model  the  fluxes,  such  as  (w/  +  w* )w„ ,  are  calculated  from  density 

functions  across  both  the  size  sections  and  vertical  cells.  In  the  model  the  size  distribution  is 
represented  by  discrete  size  sections  of  finite  width.  To  increase  the  accuracy  of  the  model  the 
mass  concentration  within  a  single  section  is  represented  by  a  density  function 


Where  v  is  the  mass  (or  volume)  of  the  particle  and  where  the  slope  <j  is  chosen  such  that  the 
number  concentration  of  particles  is  preserved  for  number  preserving  aerosol  processes  (Gelbard 
et  al.,  1998).  Integration  of  the  density  function  across  the  section  i  is,  by  definition,  the  mass 
density  m,  so  that 


«,  =- 
V 


[o-  +  l) 


m, 


CT+l 
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when  <j  =  -1 
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Where  and  are  the  lower  and  upper  boundaries  of  the  section. 

In  addition  to  the  density  function  across  the  sections,  there  is  also  a  weighting  function 
which  varies  linearly  across  the  cell  with  the  slope,  a,  determined  by  the  concentration  in 
adjacent  cells.  This  weighting  factor  is  given  by  (l  +  az)  where  z  =  0  is  the  center  of  the  cell. 
See  Gelbard  et  al.  (1998)  for  a  description  and  improved  accuracy  achieved  by  using  these 
density  functions.  The  earlier  analysis  includes  only  the  effects  of  the  gravitational  settling 
velocity.  Here  we  extend  the  analysis  to  include  upward  or  downward  fluid  motion  in  addition 
to  gravitational  settling.  We  follow  the  notation  of  the  prior  paper  in  writing  the  Stokes  fall 
velocity  as 


w®(u)  =  CD’ 


and  where  p ,  g,  and  77  are  the  density  of  the  particle,  gravitational  constant  and  kinematic 
viscosity  of  air,  respectfully*.  The  mass  removed  from  the  cell  during  a  downward  shift  is  equal 
to  the  mass  removed  between  the  bottom  of  the  cell  and  the  distance  (shift  where 

Az  is  the  width  of  the  cell  in  question.  The  change  in  mass  density  resulting  from  loss  out  the 
bottom  of  the  cell  is 


where 


c  = 


P^g 


1877 


(13) 


Am  = 


,  Az  J 


(14) 


Where  shift  is  the  vertical  shift  (as  a  function  of  size)  which  results  from  settling  and  downward 
air  movement 


shift  = 


2^ 


+co^ 


At 


Where  At  is  the  time  step.  Integration  of  Eq.  (14)  yields 


(15) 


Am  Az  \  f 

=  m< 

\-  —  {az-w^  At) 

1 

>  +  ac\ 

At  1 

[2^  J 

J 

V. 

^(az-2w^A/) 


^i+l 


4.1  N 


/+! 


caAt 

2 

( 

^  ^  J 

J 

(16) 


^  The  slip  correction  factor  has  been  neglected.  The  settling  velocity  for  particles,  which  are  small  enough  for  the 
slip  correction  to  be  important,  is  entirely  negligible  compared  to  atmospheric  motions  (turbulent  diffusion  and 
vertical  air  velocities). 
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If  =  0,  Eq.  (20)  in  Gelbard  et  al.  (1998)  is  obtained.  If  there  is  no  gravitational  settling,  c=0 
and  only  the  first  term  survives;  furthermore,  if  a  is  also  zero,  there  is  no  vertical  gradient  in  the 
mass  concentration  and  the  flux  out  of  the  cell  is  given  by  m  as  expected. 

For  a  net  upward  displacement,  i.e.  for  those  sections  where  the  upward  air  velocity 
exceeds  the  gravitational  settling,  the  limits  on  the  z-integration  in  Eq.  (14)  is  from 

( Az  I  .1^  Az 


and  the  integration  results  in  the  following  equation: 


AmAz 

At 


-ac< 


o, 


v+i 


<T+f  ^ 


cr  +  : 


+ 


caAt 

2 

^  ^  J 

(17) 


In  keeping  with  the  notation  in  the  original  paper,  a  downward  displacement  is  positive  and  an 
upward  velocity  is  negative  in  Equations  (16)  and  (17).  Equation  (17)  can  be  obtained  from  (16) 
by  (i)  changing  the  sign  on  a  and  c  and  (ii)  taking  the  absolute  value  of  . 

Equations  (16)  and  (17)  replace  the  (w^  +  w^)m„  terms  in  Equations  with  (10)  and  (1 1) 
where  it  should  be  noted  that 


AM^ii  =  (AAz)Am  =AzAm  (18) 

Without  loss  of  generality  we  can  set  A=1  m^;  i.e.,  a  column  of  unit  cross-sectional  area.  If  we 
keep  the  dimensions  of  A  we  get  aerosol  mass  density  in  units  of  mass  per  unit  volume.  In  a 
truly  1-D  situation  the  mass  density  has  units  of  mass  per  unit  length,  in  which  case  A  does  not 
enter. 


In  the  case  of  a  monodisperse  aerosol,  or  for  sections  of  infinitesimal  width,  conservation 
of  munber  of  particles  follows  directly  from  conservation  of  mass,  i.e.,  the  set  of  equations  for 

number  density  is  the  same  as  for  mass  density  and  related  by  a  single  particle  mass  for 
section  i 


N,=^  (19) 

mf 

For  the  case  where  the  mass  density  is  allowed  to  change  across  the  section,  (as  in  Eq.  [12]),  the 
number  of  particles  are  not  generally  conserved.  In  Gelbard  et  al.  (1998),  the  rigorous  method  to 
assure  conservation  of  number  as  well  as  mass  is  presented  for  the  case  of  gravitational  settling 
(no  air  motion).  However  in  implementing  particle  conservation  in  the  code,  particle 
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conservation  was  calculated  using  the  approximation  that  a  particle  which  settles  out  of  a  cell  has 
a  mass  equal  to  the  mean  particle  mass  of  the  layer  it  is  leaving.  At  the  end  of  the  time  step,  a 

new  mean  particle  mass  (new) ,  for  a  given  section,  was  calculated  from 


(new)  = 


A/r-AA^r.  AA/:, 
m^„(old)  ml^(old) 


(20) 


where  A/“"^  is  the  mass  in  the  section  at  the  start  of  the  time  step  and  the  numerator  is  the  mass 
in  cell  n  at  the  end  of  the  time  step.  The  denominator  is  the  number  of  particles  in  the  cell  at  the 
end  of  the  time  step,  (old)  is  the  mean  particle  mass,  in  the  subscripted  cell,  at  the  beginning 
of  the  time  step.  The  section  index,  i,  has  been  dropped  since  all  quantities  are  for  the  same  (dry) 
section. 


In  the  new  version  of  the  code  that  includes  vertical  air  motion,  we  have  retained  the 
approximation  that  particles  removed  from  a  cell  have  the  same  mean  particle  mass  (for  a  given 
cell)  as  the  layer  from  which  it  is  removed.  Conservation  of  number  is  now  maintained  by 
calculating  the  new  particle  mass  for  a  given  section  from 


(new)  = 


-  AAf  7'  +  AA/ ,  ±  AA/*'’'" 
-AM";'  ±AAf„"°"  AA/;„ 

m^(old)  m^±i(old) 


(21) 


where  AM""'  and  AA/,^”,  are  the  amounts  of  mass  moved  (up  or  down)  out  of  cell  n  and  into  the 
cell  n  from  adjacent  cells,  respectively.  The  direction  of  movement  is  determined  by  the  sum  of 
the  settling  velocity  and  the  vertical  air  velocity.  AA/*"'^  is  the  mass  added  to  cell  n  as 

determined  from  the  divergence  of  the  vertical  velocity  field.  Equation  (20)  is  somewhat 
schematic  and  must  be  modified  when  the  vertical  velocity  goes  to,  or  through,  zero  in  a  given 
cell. 

C.  Implementation. 

Calculations  of  organized  vertical  motions  are  carried  out  in  the  subroutine  “MOVEIT” 
(turbulent  mixing  is  done  elsewhere  in  the  code).  The  net  vertical  velocity  (settling  plus  air 
velocity)  for  a  given  section  and  cell  is  calculated  and,  depending  on  the  direction  of  the  net 
velocity,  Eq.  (16)  or  (17)  is  used  to  calculate  the  mass  removed  (AmAzAt)  and  the  mass 
remaining  in  the  cell  during  a  time  step  At.  These  masses  are  added  into  the  appropriate  cell  (and 
section)  of  a  temporary  mass  array.  The  mass  that  must  be  added/subtracted  to  each  cell  (and 
section)  as  a  result  of  convergence/divergence  of  the  fluid  field  is  calculated  from  the  last  terms 
of  Eqs.  (10)  and  (1 1).  At  the  end  of  each  time  step  (i.e.,  each  call  to  subroutine  MOVEIT)  the 
temporary  mass  array  is  used  to  re-set  the  primary  mass  array  and  to  recalculate  the  mean 
particle  mass  in  each  section  such  that  number  density  is  also  conserved. 
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For  trace  gases  the  procedure  is  much  easier  because  there  is  no  settling,  and  there  are  no 
size  sections.  There  is  however  a  density  function  across  the  cells  so  that  the  mass  flux  is  just  the 
first  term  of  Eqs.  (16)  and  (17)  with  in-flow  given  by  the  last  terms  of  Eqs.  (10)  and  (11),  which 
accounts  for  convergence/divergence. 

During  testing  it  was  found  that  a  spurious  particle  source  occurred  at  the  surface  and 
when  the  velocity  crossed  zero.  The  source  of  this  error  can  be  seen  by  writing  Eq.  (10),  for 
upward  net  motion  for  the  first  cell 

=  +wf  )/M,]+7WiW/  (22) 

Where  the  terms  for  the  cell  below  the  surface  have  been  omitted  since  nothing  can  move  up  out 
of  the  earth’s  surface.  In  the  term  that  accounts  for  convergence  there  is  no  longer  a  gravitational 
term  to  balance  the  gravitational  settling  term  in  the  cell,  as  there  is  in  all  other  cells.  This  has 
the  same  effect  as  introducing  a  large  (settling)  velocity  divergence  at  the  surface  (which  of 
course  is  not  realistic).  This  can  be  remedied  by  adding  a  gravitational  flux  term  obtained  from 
Eq.  (17)  by  setting  w  equal  to  zero  in  Eq.  (17).  This  exactly  balances  the  upward  flow  out  of  the 
first  cell  with  the  inflow;  when  there  is  no  gradient  in  mass. 


To  facilitate  one  who  is  interested  in  understanding  the  code,  we  here  provide  definitions  of 

several  key  variables  used  in  the  subroutine  MO  VEIT: 

IFLIP  is  an  integer  that  equals  1  or  -1  depending  on  whether  the  particles  in  a  given  section  and 
cell  are  moving  down  or  up,  respectively. 

NNOS  is  an  offset  index  that  is  0/2  for  downward/upward  displacement. 

TERM  is  the  mass  removed  from  the  current  section  L  and  cell  N. 

HTERM  is  the  mass  added/subtracted  (horizontally)  due  to  divergence  of  the  vertical  velocity. 

TERMOO  is  the  correction  when  the  velocity  goes  through  zero  in  the  current  cell. 

TRANQ(K,L,N)  is  the  temporary  array  for  the  mass  of  component  k  in  section  L  and  cell  N, 

where  the  array  is  offset  one  cell  (N+1)  from  that  of  the  primary  array  Q(k,l,n),  such  that 
the  mass  remaining  in  the  cell  after  the  time  step  is  placed  in  TRANQ(K,L,N+1),  the 
mass  moved  dovm  goes  in  N  and  the  mass  moved  up  goes  in  N+2. 

TRANN(L,N)  is  an  array  which  holds  the  mass  removed  vertically  from  cell  N  and  is  used  in 
calculating  the  mean  particle  mass. 

VMEANT(L,N)  is  a  temporary  array  for  VMEAN,  so  that  VMEAN(L,N)  will  not  be  overwritten 
when  moving  particles  upward. 

Analogous  variables  have  been  introduced  for  the  gas  calculations:  TRANQG(N),  HTERMG  and 
FMOVE,  the  later  being  the  mass  of  gas  removed  from  a  cell. 

In  the  program  the  sign  convention  for  the  velocity  is  positive  for  downward  velocity  and 
negative  for  upward  velocity:  and  opposite  to  that  used  in  the  preceding  equations. 

Other  definitions  remain  unchanged: 

Q(K,L,N)  is  mass  of  component  K  in  section  L  and  cell  N. 

QT(L,N)  is  the  mass  of  all  components  in  section  L  and  cell  N. 
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(The  indexing  in  the  program  is  actually  not  that  indicated  above  (K,L,N)  but  rather  a  continuous 
string:  [K+KCOMP*(L-l)+M*KCOMP*(N-l)],  where  KCOMP  and  M  are  the  number  of 
components  and  the  number  of  sections,  respectively.) 


+/-  Vertical  Velocity  (cm  ) 


0  2  4  6  8  10  12  14 

Mixing  Coef.  (m^  s“1) 

Figure  2.  Vertical  velocity  profile  (solid  line)  used  in  solutions  shown  in  Figures  3-27 
and  the  profile  of  the  turbulent  mixing  coefficient  introduced  when  mixing  is  included 

in  the  solutions 


III.  TESTING  OF  THE  MODEL 

A.  Testing  of  the  effect  of  vertical  velocity  in  isolation. 

Many  tests  of  the  effect  of  vertical  velocity  on  the  size  distribution  with  all  other  aerosol 
processes  turned  off  were  executed.  These  included:  (1)  exponentially  increasing  vertical 
velocity  with  height  to  some  height,  H,  above  which  the  velocity  remains  constant.  Both  upward 
and  downward  fluid  velocities  were  tested.  (2)  Linearly  increasing  upward  vertical  velocity 
(convergence)  to  some  height,  h,  at  which  point  the  velocity  linearly  decreased  (divergence), 
crossed  through  zero  and  became  downward,  and  at  some  greater  height,  H,  the  vertical  velocity 
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becomes  constant  in  the  downward  direction.  The  reverse  of  this  case  where  the  linearly 
increasing  velocity  is  downward  near  the  surface  was  also  checked.  (3)  Constant  upward  and 
downward  vertical  velocity  (no  convergence/divergence)  were  tested.  In  this  case  the  program 
had  to  be  modified  so  that  the  vertical  air  velocity  was  not  zero  at  the  surface  as  required  by  the 
program.  Here  we  will  show  results  only  for  the  case  (4)  where  the  vertical  velocity  increases 
linearly  with  height,  H,  after  which  it  remains  constant  with  height.  The  results  will  be  shown 
for  both  upward  and  downward  air  velocities  for  the  velocity  profile  shown  in  Figure  2  by  the 
solid  line.  For  cases  (1)  -  (3),  tests  were  limited  only  to  the  effects  of  vertical  velocity.  Many 
more  cases,  including  the  interaction  with  other  processes  were  done  for  case  (4)  some  of  which 
are  shown  below.  For  the  cases  shown  in  Figures  3-1 1 ,  turbulent  diffusion  is  turned  off  and  any 
“diffusion”  observed  is  a  result  of  numerical  diffusion.  In  later  cases  where  turbulent  diffusion  is 
included,  it  will  be  seen  to  totally  dominate  any  numerical  diffusion  which  is  introduced  by  the 
numerical  methods. 

In  all  these  tests,  it  is  useful  to  initialize  the  model  with  a  size  distribution  that  is  the  same 
at  all  altitudes  below  some  given  height  and  then  changes  abruptly  to  a  different  constant  size 
distribution  above  that  height.  The  region  of  constant  mass  with  height  gives  a  region  iii  which 
no  change  in  concentration  should  be  observed  even  if  there  is  convergence/divergence  in  the 
vertical  air  velocity  because  there  is  no  gradient  in  concentration.  The  ‘horizontal’  in/out  flow 
should  just  balance  the  difference  between  the  masses  flowing  out  the  top  and  in  the  bottom  of  a 
cell  (or  visa  versa).  The  rate  of  change  in  the  height  where  the  abrupt  change  occurs  should 
agree  with  that  predicted  by  the  vertical  velocity.  For  gases  the  motion  follows  the  air  velocity 
and  the  results  are  easier  to  interpret  than  for  particles,  the  larger  of  which,  may  have  fall 
velocities  that  are  comparable  to,  or  greater  than,  the  air  velocity. 

Downward  vertical  velocity.  The  solution  for  a  downward  vertical  velocity  profile  as  shown  in 
Figure  2  is  given  in  Figures  3-7.  All  process  are  turned  off  except  equilibration  with  the  ambient 
RH,  which  is  taken  to  be  constant  at  all  heights  at  a  value  of  70%  RH.  This  insures  that  the 
settling  velocity  of  a  given  particle  remains  constant  with  height.  The  basic  time  step  is  30 
seconds  and  data  output  every  8  hours.  Above  400  m  the  velocity  is  downward  at  1  cm  s  (no 
divergence).  Below  400  m  the  vertical  velocity  decreases  linearly  to  zero,  which  requires 
horizontal  divergence.  The  vertical  domain  of  the  model  extends  to  1550  m  (top  cell  extends 
from  1500  to  1550  m),  and  for  the  case  shown  no  material  enters  the  domain  from  the  top.  The 
thickness  of  the  vertical  cells  were  set  at  20  m  below  200  m  and  50  m  between  200  m  and  1550 
meters.  Figure  3  shows  the  evolution  of  the  concentration  profile  for  a  trace  gas  (scalar 
contaminant)  which  has  no  gravitational  settling  velocity.  The  solid  line  shows  the  initial  profile 
of  SO2.  Above  400  m  where  the  velocity  is  constant,  the  step  function  that  occurs  initially  in  the 
cell  bounded  by  600  and  650  m  and  the  discontinuity  at  the  top,  caused  by  no  new  material 
entering  the  domain,  moves  downward  at  a  constant  velocity  of  1  cm  s  .  The  discontinuity 
broadens  vertically  as  a  result  of  numerical  diffusion.  It  is  important  to  notice  that  the  profiles 
below  the  discontinuities  (and  the  diffusion  broadening)  remain  constant  with  height  -  even 
below  400  m  where  there  is  air  flow  divergence.  This  is  a  test  of  how  well  the 
convergence/divergence  part  of  the  algorithm  is  working.  While  Figure  3  is  for  a  trace  gas,  the 
same  behavior  is  obtained  for  particles  that  have  negligible  fall  velocities  (small  compared  to  0. 1 
cm  s'‘). 
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Figure  3.  Change  in  the  concentration  of  a  trace  gas  for  a  downward  vertical  velocity 
profile  shown  in  Figure  2.  The  solid  line  gives  the  initial  profile. 
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Figure  4.  Downward  movement  of  sea-salt  particles  of  dry  radius  2.4-3. 2  pm. 
Same  case  as  in  Figure  3.  Solid  line  is  initial  profile. 
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Figure  4  shows  the  downward  movement  of  the  sea-salt  component  of  particles  in  the 
section  30  where  the  dry  radius  is  2.4-3. 1  pm.  In  the  model  these  are  assigned  to  be  sea  salt 
particles  and  will  sometimes  be  referred  to  as  such.  The  downward  movement  at  the  lowest 
heights  (below  200  m)  is  greatly  enhanced  over  that  shown  in  Figure  3  as  a  result  of  gravitational 
settling.  The  solid  line  with  crosses  in  Figure  5  shows  the  initial  size  distribution  in  all  cells 
above  650  meters.  Also  shown  in  Figure  5  are  the  size  distributions  in  the  cell  at  1150-1200  m  at 
8  hrs  and  16  hours.  At  8  hours  the  size  distribution  remains  vmchanged  because  the  cells  above 
have  the  same  size  distribution  and  the  depletion  of  particles  propagating  downwards  from  the 
top  has  not  yet  reached  1200  meters.  After  about  1 1  hours  the  concentration  starts  to  drop 
precipitously  since  no  new  particles  are  being  introduced  at  the  top  of  the  domain.  The 
concentration  of  particles  after  1 6  hours  has  dropped  by  two  and  a  half  orders  of  magnitude  (it 
would  be  zero  if  it  were  not  for  numerical  diffusion). 


Radius,  (pm) 

Figure  5.  Change  in  the  size  distribution  in  the  cell  bounded  by  1150  and  1200  meters 
with  downward  velocity  profile  shown  in  Figure  2.  Solid  line  with  crosses  is  the  initial 
size  distribution  at  all  heights  above  650  meters.  No  new  particles  are  allowed  to  enter 
the  domain,  so  that  all  particles  initially  above  1150  are  below  1150  meters  afterl6 
hours. 

The  curve  denoted  by  crosses  in  Figure  6  is  the  initial  size  distribution  in  all  cells  below 
600  meters.  The  size  distribution  here  is  the  result  of  both  sulfate  and  sea-salt  aerosol,  whereas 
Figure  4  includes  only  the  sea-salt  component  which  is  dominant  at  sizes  larger  than  about  0.2 
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I^m  and  extends  initially  only  to  650  meters.^  Also  shown  is  the  evolution  of  the  size  distribution 
in  the  cell  bounded  by  400  and  450  meter.  At  16  and  24  hrs  the  initial  size  distribution  has  been 
displaced  by  the  size  distribution  which  initially  existed  above  600  meters.  In  fact,  the 
displacement  had  mostly  occurred  by  8  hours  (the  flow  displacement  in  8  hours  is  288  m).  By 
32  hours,  downward  movement  of  clean  air  from  1600  m  had  depleted  the  particle  concentration 
at  400  m  by  a  factor  of  four. 


Radius,  (inm) 

Figure  6.  Evolution  of  the  size  distribution  in  the  cell  bounded  by  400  and  450  meters 

In  the  lowest  cell  (0-20  m)  no  particles  are  removed  by  vertical  airflow  since  the  vertical 
velocity  goes  to  zero  at  the  surface.  The  removal  is  by  divergence  and  surface  deposition. 

Figure  7  shows  the  evolution  of  the  size  distribution  in  the  lowest  cell.  The  loss  of  large  particles 
is  via  gravitational  settling.  The  small  particles  moving  downward  into  the  lowest  cell  under  the 
influence  of  the  decreasing  vertical  airflow  are  just  balanced  by  horizontal  divergence  out  of  the 
cell.  The  fact  that  the  small  particle  end  of  the  size  distribution  does  not  change  with  time  is 
validation  that  the  divergence/convergence  part  of  the  program  is  working  well. 


^  The  designation  of  components  as  sulfate  and  sea-salt  aerosol  is  of  no  consequence  to  the  behavior  which  is  being 
studied  here.  However  it  is  the  way  the  components  are  specified  in  the  model.  The  aerosol  is  internally  mixed. 
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Radius,  (inm) 

Figure  7.  Evolution  of  the  size  distribution  in  the  lowest  cell  (0-20m).  The  downward 
air  velocity  is  approaching  zero  at  the  surface  and  particles  from  above  are  mostly 
removed  by  divergence  before  reaching  the  lowest  cell.  The  loss  of  larger  particles  is 
due  to  gravitational  settling. 

Upward  vertical  velocity.  The  situation  and  interpretation  is  more  complicated  for  upward  air 
velocities  because  the  net  velocity  of  the  particle  can  change  directions  depending  on  relative 
magnitudes  of  the  air  and  gravitational  settling  velocities,  and  can  even  be  in  opposite  directions 
in  adjacent  cells.  For  a  vertical  air  velocity  which  increases  linearly  with  height  there  will  be 
some  height  at  which  the  net  velocity  changes  direction.  The  numerical  solution  for  the  case 
where  the  upward  air  velocity  profile  is  given  by  Figure  2,  are  shown  in  Figmes  8  through  1 1 . 
For  particles  as  small  or  smaller  than  those  in  section  25  (0.73-0.93  mm),  the  settling  velocity  is 
so  small  that  they  move  upward  even  in  the  lowest  cell  (center  at  10  m)  as  indicated  in  the  top 
graph  of  Figure  8.  The  constancy  of  the  concentration  with  height  as  the  particles  move  upward 
indicates  that  the  (horizontal)  convergence  into  the  cells  below  400  meters  is  being  properly 
handled.  {Note:  we  will  later  set  the  air  velocity  to  zero  in  the  lowest  cell  and  discuss  the  reason 
that  it  gives  a  more  realistic  result.}  In  section  27  (1.18-1.5  pm  radius),  shown  in  the  lower  graph 
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of  Figure  8,  the  settling  velocity  in  the  lowest  cells  is  greater  than  the  upward  air  velocity  causing 
nearly  complete  depletion  of  particles  in  the  lowest  cell  in  eight  hours.  Particles  in  cells  above 
the  lowest  cells  are  moving  upward.  As  the  particle  size  increases,  the  reversal  height  moves 
upward,  as  shown  in  Figure  9,  (because  the  vertical  velocity  is  increasing  with  height).  As  the 
reversal  altitude  increases,  the  concentration  in  the  lowest  cell  is  maintained  for  an  increasing 
longer  time  by  particles  moving  down  from  above  (but  from  below  the  reversal  point)  as  is  seen 
in  Figure  9b  and  as  will  be  discussed  in  connection  Figure  10.  The  increased  settling  velocity  of 
the  larger  particles  also  slows  the  upward  movement  as  seen  in  the  figures  at  altitudes  above  650 
meters. 


Size  section  25:  0.73-.93  urn  dry  radius 


Size  section  27: 1.18-1.50  pm  dry  radius 


Sea  Salt  (10'®  kg  m'^) 


Figure  8.  Sea-salt  concentration  in  size  sections  25  and  27  at  8-hour  intervals.  Note  how 
constant  the  concentration  remains  in  the  vertical,  illustrating  that  convergence  is  being 
handled  correctly. 
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Size  section  28: 1.5-1. 9  dry  radius 


Size  section  32:  3.6-5.0  ^im  dry  radius 
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Figure  9.  Sea-salt  concentration  in  size  sections  28  and  32  at  8-hour  intervals, 
illustrating  the  increasing  height  of  the  reversal  point  as  size  of  particle  radius  increases. 

The  change  in  the  size  distribution  in  the  lowest  cell  every  eight  hours  for  48  hours  is 
shown  in  Figure  10.  The  maximums  and  minimums  developing  at  larger  sizes  can  be  explained 
by  the  behavior  shown  in  Figures  8  and  9.  The  minimum  at  about  1 .3  pm  radius  is  in  section  27 
where  the  low  concentration  results  from  the  low  reversal  point;  i.e.,  no  particle  are  available  to 
flow  downward  into  the  lowest  cell.  The  maximum  at  1.7  pm  is  in  section  28  (Figure  9a)  and 
results  because  the  reversal  point  is  now  at  a  higher  altitude  and  particles  between  the  lowest  cell 
and  the  breakpoint  are  settling  downward  into  the  lowest  cell.  At  eight  hours  all  particles  in 
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sections  with  radius  greater  than  1.5  )rm  and  smaller  than  10  (xm  are  the  same  as  at  time  zero.  Of 
course  these  sections  will  eventually  be  depleted  (as  shown  in  Figure  10)  as  the  region  between 
cell  1  and  the  reversal  point  is  depleted.  The  correct,  but  rather  strange  behavior  exhibited  in 
Figure  10  is  the  result  of  the  interplay  between  the  velocity  profile  (Figure  2),  the  dependence  of 
the  settling  velocity  on  radius,  and  the  shape  of  the  initial  concentration  and  velocity  profiles.  At 
radii  greater  than  about  5.0-6.3  pm  (section  33)  the  gravitational  settling  velocity  is  greater  than 
the  vertical  air  velocity  at  all  heights  and  the  aerosol  concentration  collapses  downward  with 
time,  as  shown  in  Figure  1 1  for  section  34. 

B.  Including  other  processes  which  effect  the  vertical  distribution. 

We  would  expect  other  processes,  which  also  cause  gradient  transport,  to  interact 
strongly  with  the  vertical  velocity  transport  discussed  above.  Other  processes  that  effect  gradient 
transport  are;  turbulent  mixing,  surface  deposition  (and/or  lower  boundary  condition),  free 
troposphere  exchange  (and/or  upper  boundary  condition),  and  a  localized  source  of  aerosols. 

Any  localized  source  of  aerosols  will  produce  a  gradient  in  concentration,  but  the  primary  one  of 
interest  here  is  the  source  of  sea-salt  aerosol  at  the  ocean  surface. 


Radius,  (urn) 

Figure  10.  Size  distribution  in  the  lowest  cell  at  8-hour  intervals 
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Size  section  34:  6.3-8. 1  urn  dry  radius 


Figure  1 1 .  Sea-salt  concentration  in  size  section  34  at  8-hour  intervals  with  upward  air 
velocity. 

Adding  turbulent  mixing.  The  vertical  broadening  (and  overshoot)  of  the  initial  step  in 
concentration,  which  occurs  between  600  and  650  m  in  Figures  8  and  9  with  increasing  time,  is 
the  result  of  numerical  diffusion.  Numerical  diffusion  is  also  apparent  in  Figure  3  for  trace 
gases.  In  a  real  atmospheric  simulation,  which  includes  turbulent  diffusion,  the  degree  of 
numerical  diffusion  seen  in  these  solutions  is  totally  obscured  by  real  turbulent  diffusion  which 
occurs  in  the  atmosphere  boundary  layer  (Fitzgerald  et  al.,  1998a) 

Figures  12  to  14  show  the  solution  when  tiubulent  diffusion  is  added,  where  the  dashed 
line  in  Figure  2  shows  the  profile  of  turbulent  mixing  coefficients  used  in  these  and  subsequent 
solutions.  This  particular  combination  of  profiles  of  turbulent  mixing  coefficients  and  aerosol 
concentrations  is  not  realistic,  since  the  drop  in  concentration  occurs  at  625  m  and  the  mixing 
coefficients  drop  off  at  about  400  m.  For  testing  the  model,  however,  it  is  better  not  to  have  the 
two  discontinuities  at  the  same  height.  Figure  12  shows  the  evolution  of  the  vertical 
concentration  profile  in  section  32,  and  should  be  compared  to  Figure  9b  when  there  is  no 
mixing.  Below  400  m  the  turbulence  is  strong  and  results  in  a  nearly  uniform  concentration. 
Above  500  m  the  mixing  drops  to  nearly  zero  and  the  upward  air  velocity  transports  the  particles 
which  were  originally  between  400  and  600  m  upward  as  in  the  case  when  there  was  no 
turbulence.  Below  about  500  m  the  combined  action  of  gravitational  settling  and  turbulence 
transports  particles  downward.  As  the  particle  concentration  tends  to  separate  at  the  reversal 
point  seen  in  Figure  9b,  the  turbulent  mixing  fills  in  the  gap  by  increasing  the  downward  flux 
from  above  and  inhibiting  the  downward  gravitational  flux  below  the  reversal  point. 

Figure  13  shows  the  case  when  the  net  (air  plus  gravitational  velocities)  is  downward  at 
all  heights.  Figure  13  should  be  compared  to  Figure  1 1  when  there  is  no  mixing.  The  mixing 
destroys  the  concentration  gradient  in  the  region  of  strong  mixing  (below  400  m).  As  the  top  of 
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the  concentration  gradient  settles  into  the  region  of  strong  mixing,  the  turbulence  acts  on  the 
gradient,  opposing  gravity,  tending  to  keep  the  particles  suspended  for  a  longer  period  of  time. 


Size  section  32:  3.6-5.0  urn  dry  radius 


Figure  12.  Size  section  32  with  up  velocity  and  turbulent  mixing.  The  turbulent  mixing 
profile  is  shown  in  Figure  2.  Figure  9b  is  the  same  section  when  there  is  no 
turbulent  mixing. 


Size  section  34:  6.3-8. 1  pm  dry  radius 


Figure  13.  Size  section  34  with  up  velocity  and  turbulent  mixing.  The  turbulent  mixing  profile 
shown  in  Figure  2.  Figure  1 1  is  for  the  same  section  when  there  is  no  mixing. 
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Figure  14  shows  the  net  effect  on  the  size  distribution  in  the  lowest  cell  over  a  48-hour 
time  period.  The  irregularities  seen  in  Figure  10  are  smoothed  out  by  the  turbulence  as  expected. 
The  one  peak  that  remains  can  be  traced  to  those  particles  that  are  nearly  suspended  by  the 
upward  air  motion  between  400  and  600  m,  above  the  region  of  mixing.  This  peak  occurs  in  size 
section  34  which  is  the  section  shown  in  FigurelS.  These  particles  fall  slowly  into  the  well- 
mixed  region  sustaining  the  particle  concentration.  Air  motion  moves  smaller  particles  upward 
and  out  of  this  region  while  larger  particles  fall  quickly  out  of  this  region. 


Radius,  (pm) 

Figure  14.  Size  distribution  at  the  surface  with  upward  velocity  and  mixing  turned  on. 

Profiles  of  mixing  coefficients  and  vertical  velocities  are  those  shown  in  Figure  2. 

Effect  of  a  sea-salt  aerosol  surface  source  and  surface  deposition.  Here  we  turned  mixing  back 
off  and  turned  on  the  surface  source  of  sea-salt  using  a  hybrid  formulation  (combination  of 
Monahan  et  al.,  1986  and  Smith  et  al.,  1993)  given  by  Hoppel  et  al.  (2002).  We  also  extended  the 
size  range  to  sizes  as  small  as  0.06  pm  even  though  there  is  no  justification  for  using  the 
formulation  below  about  0.3  pm.  Surface  deposition  was  also  turned  on.  The  surface  deposition 
was  that  given  in  Fitzgerald  et  al.  (1998a)  for  all  components  except  sea  salt.  For  sea-salt  aerosol 
we  used  only  gravitational  settling  as  shown  by  Hoppel  et  al.  (2002, 2005)  to  be  more 
appropriate  for  a  surface  source.  The  solution  uses  the  same  upward  directed  vertical  velocity 
profile  as  before  and  shown  in  Figure  2.  However,  in  the  lowest  cell  the  air  velocity  is  set  equal 
to  zero  instead  of  the  value  at  the  midpoint.  The  reason  for  this  is  that  all  particles  in  the  cell  (for 
a  given  section)  move  with  the  velocity  equal  to  the  sum  of  the  vertical  air  velocity  and  the 
settling  velocity.  For  downward  air  velocity,  this  gives  a  deposition  velocity  that  can  be  much 
too  large,  and  for  upward  air  velocity,  gives  a  deposition  velocity  that  is  either  too  small  or 
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nonexistent  if  the  upward  velocity  is  greater  than  the  settling  velocity.  Since  the  deposition 
velocity  must  indeed  go  to  zero  at  the  surface,  setting  the  air  velocity  to  zero  in  the  lowest  cell  is 
not  only  reasonable  but  totally  justifiable  by  the  fact  that  turbulent  mixing  near  the  surface  is 
much  greater  than  that  resulting  from  vertical  velocities  caused  by  large  scale  divergence.  The 
vertical  velocity  in  the  lowest  cell  is  now  set  to  zero  in  the  current  version  of  the  program. 


Size  sec.  25;  0.73  -.93  pm  dry  radius 


S-S  source  &  deposition  &  up  velocity 


Figure  15.  Vertical  profiles  of  sea-salt  concentration  in  size  section  25  when  a  sea-salt 
source  is  added  and  there  is  no  mixing. 

The  results  are  shown  for  sections  25  and  30  in  Figures  15  and  16.  In  section  25  the  net 
particle  velocity  is  upward  in  all  vertical  cells  (except  the  first  where  the  air  velocity  is  zero). 
Since  there  is  no  upward  velocity  in  the  first  cell,  the  sea-salt  aerosol  is  confined  to  the  lowest 
cell  and  the  concentration  builds  up  until  there  is  a  balance  between  the  source  and  surface 
deposition.  The  anomalous  increasing  in  the  concentration  seen  in  Figure  15  in  the  second  cell  is 
caused  by  the  following  numerical  situation:  As  stated  earlier,  the  concentration  distribution 
vertically  across  the  cell  is  assumed  to  be  linear  with  slope  determined  by  the  concentrations  in 
the  adjacent  cells.  The  unrealistic  large  build  up  in  the  first  cell  gives  a  very  large  slope  in  the 
second  cell.  Too  few  particles  are  therefore  transported  upward  out  of  the  cell  for  the  number 
moved  into  the  cell  as  calculated  from  the  velocity  divergence  using  the  midpoint  concentration. 
This  results  in  a  net  increase  in  the  concentration  in  the  cell  when  it  should  have  been  zero.  This 
is  an  unrealistic  case.  When  turbulence  mixing  is  present  the  gradient  is  small  and  nearly  linear. 
The  minimum  just  above  the  sea-salt  aerosol  layer  is  a  result  of  numerical  diffusion  that  is 
particularly  pronounced  because  of  the  very  strong  gradient  at  the  surface. 
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Size  section  30:  2.4-3. 1  urn  dry  radius 

Sea  Salt  source  (hybred)  &  deposition  &  up  velocity 


Figure  16.  Vertical  profiles  of  sea-salt  concentration  in  size  section  30  when  the  sea- 
salt  source  is  added. 

In  Figure  16,  size  section  30,  we  again  see  the  minimum  caused  by  the  reversal  in  the  net 
velocity  just  below  the  100  m  height  as  in  Figure  9b.  The  concentration  in  the  first  cell  is  limited 
by  fallout  and  is  therefore  depends  strongly  on  the  particle  radius  as  well  as  the  source  strength. 
Figure  1 7  shows  the  change  in  the  size  distribution  over  the  first  48  hours  in  the  lowest  cell.  The 
decrease  seen  at  the  smallest  radii  is  caused  by  surface  deposition.  There  is  no  sea-salt  aerosol 
source  below  about  0.06  pm  radius.  No  particle  lost  due  to  the  divergence  in  velocity  because 
the  divergence  in  vertical  velocity  is  just  matched  by  the  convergence  of  horizontal  velocity.  In 
the  region  between  0.06  and  0.5  pm  there  is  a  steady  increase  over  the  48  hour  period  because 
the  deposition  is  so  slow  that  equilibrium  (between  surface  source  and  fall  out)  is  not  achieved  in 
48  hours.  From  about  0.5  to  10  pm  equilibrium  between  the  source  and  deposition  is  achieved  in 
48  hours.  The  increasing  time  to  achieve  equilibrium  at  the  largest  radii  is  the  result  of  particles 
settling  into  the  lowest  cell  from  above. 

Add  turbulent  mixing.  Figures  1 8-20  show  the  change  that  results  when  turbulent  mixing  is 
include  where  the  profile  of  turbulent  mixing  coefficients  is  given  in  Figure  2  and  extends  to 
about  500  meters.  Comparing  Figure  18  to  15,  we  see  that  the  particles  injected  at  the  surface 
are  spread  rather  uniformly  over  the  lowest  400  m  and  reach  an  altitude  where  the  upward  air 
velocity  can  transport  them  to  higher  altitudes.  In  Figure  1 8  we  have  also  shown  the  48-hour 
solution.  Even  after  48  hours  the  concentration  is  far  from  the  steady-state  value  (between 
upward  gradient  mixing  and  gravitational  fallout).  The  time  to  reach  equilibrium  is  inversely 
proportional  to  the  fall  velocity  and  for  small  particles  equilibrium  is  never  reached  in  the 
atmosphere  (Hoppel,  2002). 
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Cell  1:  Upward  vel.  &  Sea-salt  source 
&  surface  deposition 


Radius,  (|im) 

Figure  17.  Evolution  of  the  size  distribution  in  the  lowest  cell  over  48  hours.  Processes 
include  a  sea-salt  surface  source  and  surface  deposition  (and  no  mixing) 

The  newly  injected  particles  in  section  30  (2.4-3. 1  pm  dry  radius),  are  no  longer  confined 
to  the  lowest  cell  as  shown  in  Figure  1 6,  but  are  now  mixed  throughout  the  first  400  m  as  seen  in 
Figure  19.  Many  of  the  particles,  initially  in  the  upper  part  of  the  400  m  layer,  are  mixed  down 
to  the  surface  where  they  are  subject  to  removal  by  gravitational  settling.  (Figures  19  and  16 
have  different  vertical  scales)  The  equilibrium  surface  concentration  for  Section  30  is  about  half 
the  original  concentration. 

Figure  20  shows  the  evolution  of  the  size  distribution  in  the  lowest  cell  over  a  48  hour 
period  when  turbulence  is  presence.  The  large  depletion  of  particles  at  small  radii  due  to  surface 
deposition  seen  in  Figure  17  has  been  spread  over  400  meters  greatly  reducing  the  depletion  in 
cell  1 .  Likewise  the  large  increase  at  intermediate  sizes  due  to  sea-salt  injection  into  the  lowest 
cell  is  spread  to  a  greater  vertical  extent  reducing  the  concentration  in  the  lowest  cell.  Particles 
larger  than  about  1  pm  appear  to  reach  an  equilibrium  concentration  that  is  less  than  the  initial 
concentration  within  the  48  hour  time  period.  The  concentration  of  particles  smaller  than  1pm 
(and  larger  than  the  sea-salt  source  cutoff)  keeps  increasing  because  of  the  extremely  long  time 
to  reach  equilibrium,  which  is  the  result  of  the  very  small  settling  velocity.  The  time  constant  for 
achieving  equilibrium  at  radii  of  about  1  to  3  pm  inferred  from  Figure  20  is  much  faster  than  the 
actual  time  constant  and  that  given  in  Hoppel  et  al.  (2002).  In  this  particular  case,  we 
fortuitously  started  with  concentrations  very  near  their  equilibrium  values  in  this  size  range  as 
indicated  by  the  concentration  at  about  1  pm  radius  in  Figure  20  which  does  not  change  with 
time.  The  size  distributions  (not  shown)  in  all  cells  below  500  m  and  at  radii  smaller  than  5  pm 
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look  almost  identical  indicating  that  the  atmosphere  is  well  mixed  to  that  height.  For  particles 
larger  than  about  5  pm,  the  concentration  drops  off  more  rapidly  with  increasing  size  as  the 
height  increases,  reflecting  the  fact  that  gravitational  settling  is  sufficiently  strong  to  establish  a 
vertical  gradient  in  the  concentration  of  large  particles. 


Size  sec.  25:  0.73  -.93  pm  dry  radius 
S“S  source  &  deposition  &  up  velocity&  mix 


Figure  18.  Change  in  the  vertical  profile  in  section  25  when  turbulence  is  added. 
Compare  with  Figure  15. 

Size  section  30: 2.4-3. 1  pm  dry  radius 

S-S  source  &  deposition  &  up  velocity  &  mix 


Figure  19.  Change  in  the  vertical  profile  in  section  30  when  turbulence  is  added. 
Compare  with  Figure  16. 
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Cell  1 :  upward  vel.  &  Sea-salt  source 


Radius,  (^m) 

Figure  20.  Evolution  of  the  size  distribution  when  turbulence  is  added.  Compare  with 
Figure  17. 


Radius,  (pm) 

Figure  21.  Evolution  of  the  size  distribution  at  1150-1200  meters  when  turbulence  is 
added 
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Figure  21  shows  the  evolution  of  the  size  distribution  in  the  cell  at  1150-1200  meters.  At 
this  height  turbulent  mixing  is  negligible  and  the  size  distribution  remains  constant  for  8  hours 
because  the  air  between  600  and  1 150  m  has  the  same  initial  size  distribution.  By  16  hours  the 
air  from  the  lower  layers  has  reached  this  cell  and  the  size  distribution  looks  much  like  that  in  the 
mixed  layer  except  for  the  gravitational  filtering  out  of  the  largest  particles. 

Downward  velocity.  The  prior  solutions  were  for  an  upward  directed  velocity.  Figure  22  shows 
the  solution  for  the  size  distribution  in  the  lowest  cell  for  the  case  of  a  downward  velocity  with 
profile  shown  in  Figure  2.  The  processes  included  are  the  same  as  the  prior  case  shown  for 
upward  velocity  in  Figure  20.  The  interpretation  of  Figure  22  is  quite  instructive.  For  the  first  8 
hours  there  is  no  change  except  that  caused  by  gravitational  settling  at  radii  greater  than  about  5 
[xm.  After  8  hours  subsidence  brings  air  from  above  down  into  the  mixed  layer,  and  the  size 
distribution  at  radii  below  about  0.1  pm  approaches  that  of  the  FT  at  16, 24  and  32  hours.  (The 
initial  size  distribution  above  the  mixed  layer  is  shown  in  Figure  5  and  again  in  Figure  23  as 
square  symbols).  At  40  and  48  hours  subsidence  has  brought  particles  free  air  fi-om  above  the 
domain  into  the  top  of  the  mixed  layer,  which  causes  the  dilution  seen  at  radii  smaller  than  about 
0.1  pm.  At  larger  sizes  we  see  continuous  dilution  after  8  hours,  since  there  are  no  large 
particles  in  the  upper  air.  However,  the  continuous  injection  of  sea-salt  particles  at  the  siuface 
causes  equilibrium  to  be  achieved  in  the  mixed  region.  As  seen  in  Figure  22,  the  time  required  to 
reach  equilibrium  decreases  with  particle  radius.  At  1 0  pm  radius  the  equilibrium  is  reached  in  8 
hours,  but  requires  48  hours  at  1  pm. 

If  we  now  let  the  air  above  the  domain  (which  we  will  here  refer  to  as  FT  air)  have  the 
same  size  distribution  as  the  upper  air  and  allow  exchange,  as  determined  by  an  exchange 
velocity.  Vex ,  we  get  the  evolution  of  the  size  distribution  shown  in  Figure  23.  One  might  be 
surprised  that  the  magnitude  of  the  size  distribution  at  small  radii  after  48  hours  is  only  about 
half  that  of  the  FT  size  distribution.  In  the  case  shown  in  Figure  23  we  have  used  an  exchange 
velocity  of  1  cm  s'*,  which  is  the  same  as  the  downward  velocity  in  the  top  cell.  In  steady  state 
the  flux  into  the  top  of  the  domain  must  equal  the  flux  out  of  the  bottom  of  the  top  cell 


GnccIi)  ~^NceuQNcell 


GnccII 


V  +V 

'  ex  ^  '  Ncell 


QFT 


where  OnccII  and  QFT  are  the  mass  concentrations  in  the  top  cell  and  above  the  domain 
respectively.  Figure  24  shows  the  evolution  of  the  profile  of  mass  in  a  small  size  section  where 
gravitational  settling  is  negligible.  The  equilibrium  value  is  half  the  FT  value  in  keeping  with 
the  above  equation.  The  mixing  below  about  400  m  keeps  the  gradient  in  concentration  below 
400  m  small.  To  get  the  concentration  to  remain  at  its  initial  value  would  require  that  Vex» 
Vnccii-  (For  larger  particles  VNceii  includes  gravitational  settling.) 
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Cell  1:  Downward  vel.  &  Sea-salt  source 


Radius,  (^m) 

Figure  22.  Evolution  of  the  size  distribution  over  48  hours  when  the  direction  of  the 
vertical  velocity  profile  is  downward.  The  sea-salt  aerosol  source,  surface  deposition 
and  mixing  are  all  on. 


Cell  1 :  Downward  vel.  &  Sea-salt  source 


Radius,  (pm) 

Figure  23.  Same  as  in  Figure  22  but  with  free  troposphere  exchange  also  turned  on,  where  the 
free  troposphere  size  distribution  is  shown  the  blue  squares. 
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Size  section  7:  .019-.025  |am  dry  radius 

S-S  source  &  deposition  &  dn  velocity  &  mix  &  FT  Ex. 


Figure  24.  Time  evolution  of  Section  7  when  FT  exchange  is  on  and  the  exchange 
velocity  is  equal  to  the  downward  air  velocity  in  the  top  cell. 


The  model  allows  the  mechanism  of  exchange  with  the  FT  to  be  either  through  an 
exchange  velocity  as  illustrated  above  or  by  a  downward  (subsidence)  velocity  into  the  top  cell 
(or,  as  illustrated  earlier,  an  upward  velocity  out  of  the  domain).  An  exchange  velocity  is 
sometimes  provided  from  experimental  data  (usually  derived  from  moisture  or  ozone  exchange). 

If  instead  of  using  the  exchange  velocity  (gradient  transfer)  to  interact  with  the  FT,  we 
allow  an  downward  airflow  to  bring  mass  in  from  above  we  get  the  solution  shown  in  Figures 
25-27.  The  downward  velocity  into  the  top  cell  from  above  is  vertical  air  velocity  plus  the 
settling  velocity  for  particles  and  just  the  air  velocity  for  trace  gases.  For  the  solutions  shown  we 
have  set  the  mass  concentrations  above  the  domain  to  be  the  values  in  the  top  cell  of  the  domain. 
In  this  case  there  should  be  no  temporal  change  in  the  concentration  in  the  top  cell.  Hence,  in 
Figure  25,  subsidence  causes  no  change  in  concentration  in  the  upper  cells.  In  the  mixed  layer 
the  concentration  approaches  the  FT  concentration  as  subsiding  air  displaces  the  original  air. 
Figure  25  shows  that  the  small  radius  end  of  the  size  distribution  approaches  the  FT  distribution 
and  the  large  salt  particles  mix  upward,  and  strike  a  balance  between  the  source  and  surface 
deposition  of  particles  as  before.  The  profile  of  DMS,  shown  in  Figure  27  is  interesting.  There 
is  no  DMS  initially  below  600  m,  and  a  low  concentration  above.  Subsiding  air  brings  DMS 
downward  and  the  concentration  decreases  somewhat  because  of  oxidation  by  OH.  DMS  is 
mixed  upward  from  the  surface,  but  DMS  from  the  surface  source  cannot  penetrate  above  500 
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meters,  the  extent  of  the  mixed  layer.  The  DMS  concentration  in  the  mixed  layer  is,  under  these 
circumstances  determined  by  the  source  at  the  surface,  surface  deposition,  dilution  by  subsiding 
air  from  above,  and  loss  resulting  from  OH  oxidation. 


Size  section  7;  .01 9-. 025  pm  dry  radius 
S-S  source  &  deposition  &  dn  velocity  &  mix 


Figure  25.  Same  as  in  Figure  24  except  FT  exchange  is  turned  off  and  material  flows  in 
from  above  with  vertical  downward  velocity  in  the  top  cell. 
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dN/d[log  r],  (cm 


Cell  1 :  Downward  vel.  &  Sea-salt  source 
&  surf.  dep.  &  turb.  mix  &  F.T. Exchange 
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Figure  26.  Evolution  of  the  size  distribution  in  the  lowest  cell  when  there  is  subsidence  into  the 
top  from  above.  At  small  radii  the  MBL  size  distribution  is  determined  by  the  subsiding  air  and 
at  large  sizes  by  the  surface  source  of  sea-salt  aerosol. 
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Dn  velocity,  DMS  source  =  5.74X10'^^ 


Figure  27.  A  small  concentration  of  DMS  is  supplied  from  above  by  subsidence  where 
the  decrease  in  the  subsiding  air  is  due  to  oxidation  ay  OH.  In  the  mixed  layer  DMS  is 
mixed  upward  from  the  surface  source. 


IV.  CONCLUSIONS 

The  above  tests,  plus  ones  not  shown,  demonstrate  that  the  addition  of  large-scale 
atmospheric  motions  (subsidence  and  lifting)  can  be  added  to  the  NRL  aerosol  and  trace  gas 
model  (MARBLES).  The  large-scale  convergence  and  divergences  that  cause  the  vertical 
motion  must  be  generated  externally.  The  current  source  of  the  vertical  velocities  is  the  Navy’s 
3-D  mesoscale  model  COAMPS  (Coupled  Ocean  Atmosphere  Meteorological  Prediction 
System).  The  aerosol  model  is  run  along  an  air-mass  trajectory  generated  from  the  output  of 
CO  AMPS  which  also  supplies  all  the  vertical  profiles  of  meteorological  data  required  by  the 
aerosol  model  along  the  trajectory  (i.e.;  turbulent  mixing  coefficients,  large-scale  vertical 
velocities,  temperature,  relative  humidity,  etc.).  The  1-D  aerosol  model  can  be  visualized  as  a 
Lagrangian  column  moving  with  a  mean  horizontal  wind  speed  of  the  MBL.  The  1-D  limitation 
requires  horizontal  homogeneity,  but  the  vertical  stratification  moves  upward  or  downward,  or  is 
compressed  or  expanded  by  large-scale  subsidence  and  lifting.  This  (1-D)  limitation  arises 
because  the  air  that  enters  or  leaves  the  column,  at  any  designated  height,  must  have  the  same 
aerosol  load  both  inside  and  outside  the  column.  Within  any  cell  (horizontal  layer),  all  the 
aerosol  source,  sinks  and  transformation  processes,  which  are  responsible  for  the  evolution  of  the 
size  distribution,  are  taking  place. 
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The  tests  shown  above  were  designed  to  test  the  algorithms  and  therefore  do  not  represent 
realistic  atmospheric  conditions.  Nevertheless,  working  through  the  examples  will  not  only  help 
one  to  understand  how  the  numerical  simulations  work,  but  will  also  contribute  to  a  better 
understanding  of  how  atmospheric  dynamics  affect  the  vertical  variations  of  aerosols. 

For  an  ideal  Lagrangian  column  there  is  no  exchange  through  the  walls  of  the  column. 
The  modifications  introduced  by  the  work  reported  here  allow  transport  through  the  walls  due  to 
convergence  and  divergence  of  the  local  flow  field,  with  the  restriction,  that  the  air  entering  the 
column  at  a  given  height  has  the  same  aerosol  load  inside  and  outside  the  eolumn  (this  is  the  1  -D 
or  horizontally-homogeneous  limitation).  It  is  clear  that  in  regions  where  there  are  localized 
sources  and  sinks,  the  1-D  assumption  will  not  be  valid.  It  is  useful  to  further  examine  the  1-D 
limitation  and  define  under  what  conditions  the  1-D  assumption  holds. 

First  consider  the  case  of  large-scale  subsidence  causing  horizontal  divergence  in  the 
MBL  flow  field.  In  this  case  the  airflow  will  be  out  the  sides  of  the  column  and  horizontal 
inhomogeneities  outside  the  column  will  not  affect  conditions  within  the  column.  The  subsiding 
air  will  propagate  downward  through  the  top  of  the  column  and  be  mixed  with  the  preexisting  air 
by  turbulent  diffusion.  In  this  case,  the  1-D  limitation  does  not  appear  to  be  significant,  but  the 
accuracy  with  which  we  can  specify  the  flux  into  the  top  of  the  column  is  critical.  The  vertical 
variation  of  large  sea-salt  particles  (radius  greater  than  about  5  pm)  will  occur  more  rapidly  than 
dilution  by  mixing  so  that  the  subsidence  will  have  little  effect  on  large  sea-salt  particles. 

Smaller  particles,  whose  lifetime  is  much  longer  (determined  primarily  by  precipitation 
scavenging),  will  be  well  mixed  vertically  and  the  concentration  will  slowly  change  to  reflect 
that  in  the  air  above  the  MBL. 

For  the  case  of  convergence  within  the  MBL,  with  upward  velocity  in  the  column, 
horizontal-inhomogeneities  outside  the  column  can  be  advected  into  the  column.  Except  for 
large  sea-salt  particles  where  the  establishment  of  the  gravitationally  induced  vertical  gradient  is 
more  rapid  than  the  convergence,  the  potential  for  change  within  the  column  caused  by 
convergence  is  real.  To  further  evaluate  the  effect  of  horizontal  inhomogeneities  during 
convergence  consider  the  following  situation. 

(1)  Consider  a  Lagrangian  column  that  is  being  advected  along  an  air  mass  trajectory.  In 
principal,  a  Lagrangian  column  is  defined  such  that  the  mean  horizontal  flow  around  the 
periphery  is  zero  (except  for  the  small  amount  due  to  local  convergence).  Since  the  model  grid 
spacing  is  the  order  of  a  kilometer,  we  take  the  radius  ri  of  the  column  as  1  km.  The  Eulerian 
wind  field  of  the  mesoscale  model  includes  the  effects  of  the  large-scale  convergence.  The  local 
convergence  in  the  Lagrangian  flow  field  around  the  column  is  what  gives  rise  to  the  vertical 
velocity  in  the  column. 

(2)  Further  consider  a  cylinder  of  radius  ro  and  height  Hmw  as  shown  in  Figure  A1  coaxial  with 
the  Lagrangian  column.  The  region  inside  ro  is  assumed  to  be  horizontally  homogeneous.  The 
region  outside  to  is  horizontally  inhomogeneous  and  we  wish  to  calculate  the  time  required  for 
the  inhomogeneous  region  at  to  to  reach  the  column  rj  for  an  upward  velocity  Uup  which  is 
uniform  over  the  upper  surface  at  Hmw  >  where  Hmw  is  the  height  of  the  MBL.  The  mean  radial 
velocity  is  given  by 
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The  time  to  travel  from  ro  to  r,  is  then 


(A2) 


If  the  horizontally  homogeneous  region  has  a  radius  of  10  km  (r,,  =  70  km)  and  if  the  MBL  is 
assumed  to  have  a  height  of  1000  m  with  uniform  lifting  of  1  cm  s*',  then  the  time  for  the 
inhomogeneities  to  reach  the  column  is  about  5  days.  This  is  about  the  same  time  as  a  mean 
lifetime  for  an  aerosol  particle  at  the  (accumulation  mode)  radius  where  the  maximum  lifetime  in 
the  MBL  occurs  (Hoppel  et  al.  [1990],  Figure  24).  If  the  time  for  the  local  convergence  to 
advect  the  aerosol  from  the  region  of  the  inhomogeneity  to  the  column  is  greater  than  the  mean 
lifetime  of  the  aerosol,  then  other  dynamic  processes  already  accounted  for  in  the  column  can 
resolve  the  evolution  of  the  size  distribution.  We  know  from  measurements  that  horizontal 
homogeneity  can  extend  over  large  areas  of  the  ocean  with  more  rapid  changes  occurring  in 
frontal  areas  and  near  coasts  (see,  for  example,  Hoppel  et  al.  [1990]).  While  the  above  analysis 
is  very  rough  it  does  give  us  a  clue  as  to  when  we  would  expect  the  Lagrangian  column  model, 
as  modified  here,  to  give  representative  results.  With  regard  to  the  1-D  limitation,  the  model  is 
expected  to  be  more  representative  when  the  column  is  advected  through  regions  of  subsidence 
(high  pressure  regions)  than  regions  of  convergence.  In  regions  of  convergence  (vertical  lifting) 
horizontal  gradients  of  aerosol  concentrations  must  be  weak  over  regions  the  order  of  10  km 
around  the  advecting  column  for  the  1-D  approximation  to  be  valid. 

We  have  not  mentioned  horizontal  mixing  due  to  turbulence  as  a  mechanism  for 
transporting  horizontal  gradients.  Estimates  indicate  that  horizontal  mixing  in  the  Lagrangian 
frame  over  distances  the  order  of  10  km  is  not  as  severe  as  the  problem  discussed  above. 
Turbulent  mixing  has  a  tendency  to  destroy  horizontal  gradients. 


Hmbl 


Figure  28.  r;  is  the  radius  of  the  Lagrangian  column,  ro  is  the  radius  whieh  divides  the 
homogeneous  from  the  inhomogeneous. 
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APPENDIX:  Including  the  density  stratiflcation  of  the  atmosphere. 

The  analysis  given  in  this  report  can  and  should  be  expanded  to  include  the  density 
stratification  of  the  atmosphere.  Under  typical  atmospheric  conditions  the  density  decreases  by 
about  1 0%  in  the  first  1 000  m  and  1 8%  in  the  lowest  2000  meters. 

Let  p  be  the  air  density,  then  conservation  of  mass  requires 

^  =  (Al) 

dt 

For  the  purposes  of  calculating  the  density  stratification,  the  atmosphere  can  be  considered  to  be 
in  static  equilibrium  dp  j 5/  =  0 . 

Following  the  same  procedure  that  led  to  Eq.  (7),  we  obtain  Eq.  (A2)  from  (Al)  for  the 
case  of  a  density-stratified  atmosphere. 


-  ppv^A  +  pp^k^  -  •dS  =  Q 

hor 

Combining  Eq.  (A2)  and  Eq.  (5),  we  find  the  change  of  the  mass  in  a  cell  to  be 


(A2) 


A  dt 


-{wl  +wl 


(A3) 


where  (yo)  is  the  mean  density  in  the  cell  which  can  be  taken  as  the  density  at  the  center  of  the 
cell.  Equation  (A3)  replaces  Eq.  (8),  which  was  for  an  incompressible  atmosphere.  Values  for 
Pi,  P^,  and  (p) ,  for  a  given  cell  can  be  obtained  from  the  meteorological  model,  sounding,  or 

standard  atmosphere  as  the  application  warrants.  Equation  (A3)  predicts  a  change  in  aerosol 
(sectional)  mass  even  if  there  is  no  velocity  gradient  and  no  aerosol  mass  gradient 


A  dt 


(a -a) 


(A4) 


where  and  m  are  constant  with  height,  and  5  is  the  fi'actional  change  in  density  across  the 
cell. 
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It  may  seem  contradictory  to  have  a  change  in  the  aerosol  mass  in  a  cell  when  the  aerosol 
mass  flux  is  the  same  at  the  top  and  bottom  of  the  cell.  However  the  difference  in  m  mass  flux 
at  the  top  and  the  bottom  requires  an  inflow  of  air  (horizontally)  given  by  Eq.  (A4).  If  we  want 
the  mass  in  the  cell  to  remain  constant  when  the  vertical  velocity  is  constant  with  height,  Eq. 

(A3)  requires  that  the  fractional  change  in  the  aerosol  mass  density  be  equal  (and  opposite)  to  the 
fractional  change  in  the  air  mass  density. 

{”’)  (p) 

where  we  have  assumed  that  the  gravitational  settling  velocity  does  not  change  with  height  (the 
case  with  no  change  in  relative  humidity  with  height). 
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